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ABSTRACT 


The  internal  heat  transfer  problem  for  a  typical 
power-transistor  structure  has  been  solved  analyti¬ 
cally.  The  relations  among  current  distribution, 
heat  generation,  and  temperature  distribution  have 
been  derived.  Usage  of  the  resulting  equations  is 
illustrated  by  application  to  the  most  elementary 


problem, 


namely,  uniform  heat  generation  under 


the  emitter. 
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I.  INTRODUCTION 


One  limiting  factor  in  the  operation  of  a  transistor  is  the  internal  heating 
due  to  the  flow  of  current.  Not  only  does  heating  affect  the  characteristics  of 
the  transistor,  but  in  excessive  amounts  it  may  also  cause  certain  junction 
areas  to  become  nearly  intrinsic.  When  this  happens,  the  transistor  action 
ceases  altogether.  In  power  transistors,  which  are  intended  to  carry  large 
current  loads,  the  internal  heating  is  an  important  design  consideration.  In 
this  paper,  the  heat  transfer  characteristics  of  a  typical  power  transistor  are 
examined  theoretically. 


II.  THEORETICAL  ANALYSIS 


A.  BASIC  SOLUTION 

A  diagram  of  one  type  of  power  transistor  structure  is  shown  in  Fig.  1. 
On  the  top  of  the  transistor  is  a  parallel  series  of  long  emitters.  The  emitters 
and  the  base  region  are  actually  much  thinner  than  shown  in  the  side  section. 
The  collector  is  securely  attached  to  a  heat  sink.  Most  of  the  heat  generated 
in  the  structure  is  developed  in  the  high- resistivity  region  near  the  collector- 
base  junction.  Virtually  all  of  the  heat  generated  in  the  transistor  is  removed 
by  the  heat  sink.  The  problem  is  to  derive  an  expression  for  the  temperature 
distribution  inside  the  structure  in  order  to  locate  any  hot  spots  in  the  junction 
regions. 

It  was  expected  that,  because  of  end  effects,  the  hottest  portion  of  each 
emitter  was  half-way  down  its  length  (topviewin  Fig.  1);  hence,  we  solved  for 
this  area  only.  A  cross-section  of  a  single  emitter  was  examined,  using  two 
of  several  possible  techniques.  One  of  the  techniques  used  was  to  solve  for  a 
single  heat  source  on  one  side  of  infinite  parallel  planes  (thus  neglecting  side 
effects).  In  this  case,  the  total  solution  for  several  emitters  would  be  obtained 
merely  by  superposing  the  individual  solutions.  Unfortunately,  however,  we 
were  unable  to  obtain  the  desired  result  by  using  this  technique. 


Fig.  1.  Structure  of  Typical  Power  Transistor 

Collector  is  actually  very  much 
thicker  than  base  plus  emitter. 


EMITTER  MIDPOINT 

CENTER  BETWEEN 


y  *  b,  T  •  T, 


Fig.  2.  Geometry  and  Coordinate  System  for  Analysis  of 
Heat  Transfer  in  Power  Transistor  of  Fig.  1 
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As  an  alternate  approach,  we  considered  the  central  emitter  in  a  large 
collection.  The  central  emitter  was  expected  to  be  the  warmest;  hence,  this 
was  the  one  of  most  interest.  By  considering  the  central  emitter,  we  again 
neglected  side  effects;  however,  with  this  technique,  the  final  result  should 
still  be  quite  good  even  if  the  total  number  of  emitters  is  as  small  as  five. 

A  model  of  the  problem  under  consideration  is  shown  in  Fig.  2.  Because 
of  symmetry,  no  heat  flows  across  the  dividing  plane  half-way  between  the 
emitters  and  no  heat  flows  from  one  side  of  an  emitter  to  another.  Hence, 
these  two  boundaries  are  considered  to  be  insulated.  The  heat  sink, on  the 
bottom  of  the  rectangle  at  y  =  b,  is  at  temperature  Tg.  Half  of  the  collector- 
base  junction  beneath  an  emitter  occupies  the  uppe  r  left  portion  of  the  rectangle 
at  y  =  0. 

The  assumptions  made  in  this  analysis  are  summarized  as  follows: 


.  * 

is 


1)  homogeneous  and  isotropic  thermal  conductivity 

2)  steady  state  conditions 

3)  infinite  number  of  infinitely  long  emitters,  i.  e.  ,  end 
effects  are  neglected 

4)  heat  is  generated  only  along  the  collector-base  junction, 
y  =  0.  (This  primarily  requires  very  thin  emitter  and 
base  regions.) 

5)  temperature  at  y  =  b  is  uniform  at  Tg 


With  these  assumptions,  the  differential  equation  governing  this  process 


(1) 


Definitions  of  symbols  are  given  in  the  Nomenclature  at  the  end  of  this 
report. 
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and  the  boundary  conditions  are: 


8T 

at  x  =  0,  =  0 

t  8T  . 

at  x  “  a*  w =  0 


at  y =  0i  I7 =  ■ 


at  y  =  b,  T  =  T 


It  may  easily  be  shown  that  the  following  equation  satisfies  the  above  conditions, 
and  is  thus  the  general  solution  to  our  problem. 


00 


T  • T.  ■  T<b  -  t2*»(sr) 


mrx  .  ,  -  y  /  aim  , 

cos -  sinh  ■■■■■  —  sech  ■ 

a  a  a 


(bJLx)n3T  bnir 


(2) 


n=l 


The  constants  an  are  found  from  the  boundary  condition  at  y  =  0, 


3T 

=  - 


00 


n-1 


mrx  gM  . 
cos  =  -  f(x) 


This  forms  a  Fourier  series  and  hence 


a 


n 


nirx  , 

cos - dx 

a 


(n  =  0,  1,  2,  •  •  •) 


(3) 


(4) 


This  solution  could  also  be  applied  to  the  whole  device,  as  shown  in 
Fig.  1,  rather  than  to  a  single  emitter.  This  is  true  because  the  device  is 
generally  much  wider  than  it  is  thick  so  that  the  edges  may  be  considered  to  be 
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insulated.  The  only  alteration  in  the  solution  would  be  in  the  definition  of  a, 
which  would  then  become  the  total  device  width.  The  heat  generation  f(x) 
would  in  this  case  be  periodic  with  the  emitter  spacing.  A  treatment  of  this 
type  would  also  take  into  account  the  side  effects  of  a  finite  number  of  emitters. 

B.  RELATION  TO  EMITTER  TEMPERATURE 

From  an  experimental  standpoint,  temperature  is  much  more  easily 
measured  than  heat  flux.  Therefore,  it  is  of  interest  to  have  a  solution 
for  the  problem  with  the  temperature  distribution  g(x),  rather  than  the  heat 
flux,  specified  on  the  collector-base  junction.  From  Eq.  (2),  we  find  at 
y  =  0  that 


oo 

T  n  -  T  =  ^  b  +  V*  a  (— ) 
y=0  s  2  '  ^  n\nir/ 


nirx  .  .  bmr 

cos - tanh  - 

a  a 


n=  1 


(5) 


We  now  define 


b  =  a  tanh  and  b  =  a  b 

n  n\niT/  a  o  o 

Substituting  these  into  Eq.  (5),  we  obtain 


(6) 


oo 


-  T  = 
s 


o  ,  X  s ,  nirx 

— —  +  7  b  cos - 

2  n  a 


g(x) 


n=  1 


(7) 


Hence, 


f 


b 

n 


a 

2 


g(x)  cos 


mrx 

a 


dx 


(8) 
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and 


T  -  T.  ■  t(>  ■  S)  *  E bn  “■  ¥  •“»  “•“*  ¥  ,9> 

n=l 

Suppose  that  g(x)  (the  temperature  at  y  -  0)  were  measured  in  some 
way.  From  this,  f(x)  may  be  found  by  using  Eqs.  (3),  (6),  and  (8),  to  obtain 

/•a  oo 

*(x)  =  35  j  g(x)  631  +  ^(— )  cos  coth  I  g(x)  cos  dx  (10) 

Jo  n=l  Jo 

If  it  is  assumed  that  the  surface  of  the  transistor  is  at  g(x),  then  the 
temperature  could  conceivably  be  measured  by  using  a  series  of  materials 
with  progressively  increasing  melting  points.  Each  material  would  be  coated 
in  turn  onto  the  device  and  the  position  of  melting  observed. 

C.  RELATION  BETWEEN  CURRENT  DISTRIBUTION  AND 
HEAT  GENERATION 

Once  f(x)  is  obtained,  perhaps  in  the  manner  indicated  in  the  previous 
section,  it  should  be  possible  to  say  something  about  the  distribution  of  current 
under  the  emitter.  Conversely,  if  the  current  distribution  is  known,  it  should 
then  be  possible  to  arrive  at  f(x). 

It  is  assumed  that  the  heat  generation  at  any  point  under  the  emitter  is 
proportional  to  the  current  density  j(x)  at  that  point,  or  that 

j(x)  =  ^  (11) 

where  K  is  a  constant  that  will  be  determined.  This  is  of  course  only 
a  first-order  approximation;  K  would  also  be  a  weak  function  of  the 
transistor  operating  conditions. 
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If  the  total  current  per  unit  emitter  length  is  i,  then 


Hence, 


i  =  2  j  j(x)  dx  =  2  j 


K  =  2  J  ^ 


■  i  \  i  f(x) 

=  , 

2  I  f(x)  dx 

»  r\ 


In  practice,  however,  a  heavy  metallic  contact  layer  lies  over  the  emitter 
and  tends  to  even  out  the  temperature  of  the  emitter  and  invalidate  Eq.  (11). 

III.  EXAMPLES 

A.  SEMI-INFINITE  ONE -DIMENSIONAL  TRANSISTOR 

Fletcher  has  given  an  approximate  solution  for  the  current  distribution 
in  an  idealized  semi-infinite  one -dimensional  transistor.  The  solution  is 
of  the  form 


JW  -  jjl  +  AxJ' 


*N.H,  Fletcher,  "Some  Aspects  of  the  Design  of  Power  Transistors," 
Proc.  IRE  43,  551-59(1955). 
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where  A  is  a  constant  and  x@  is  the  distance  in  from  the  edge  of  the  emitter, 
w  -  x.  This  expression  tells  us  that  the  current  is  concentrated  near  the 
edges  of  the  emitters.  Unfortunately,  when  this  equation  is  substituted  into 
Eq,  (13)  and  this  in  turn  substituted  into  Eq.  (4),  the  result  cannot  be  inte¬ 
grated  analytically. 

Thus,  we  are  reminded  that,  if  analytical  results  are  desired,  f(x) 

must  be  of  such  a  form  that  f(x)  cos  mrx/a  can  be  integrated.  Furthermore, 

because  the  current  is  known  to  concentrate  near  the  edges  of  the  emitters, 

f(x)  must  be  constant  or  continually  increasing  with  x.  Some  functions  that 

3x 

fulfill  these  two  conditions  ares  A,  Ax®(s=l,  2,  .  ,  .),  A  sin  Bx,  Ae°  ,  and 
any  sums  thereof.  Because  of  its  simplicity,  only  the  application  of  the  first 
function  will  be  considered  in  detail  here. 

B.  UNIFORM  HEAT  GENERATION 

We  consider  now  the  case  of  uniform  heat  generation  beneath  an  emitter 
of  width  2w.  Let  q  be  the  rate  of  heat  generation  per  unit  length  of  emitter 
(e.g.  ,  watt/cm).  At  the  collector-base  junction  y  =  0,  we  then  have 


uaT  q 
_k57  "  Iw 

or 

. aT  -  q 

”  FjT  2kw 

Hence 

f(x)  =  ife 

for 

0  <  x  <  w 

f(x)  =  0 

for 

w  <  x  <  a 

By  substituting  this  into  Eq.  (4),  we  obtain 


(16) 


and 


Substitution  of  (16)  and  (17)  into  Eq.  (2)  then  yields  the  desired  solution  for 
uniform  heat  generation  under  the  emitters, 


CD 


In  the  actual  operation  of  the  transistor,  we  would  like  to  know  or  to 
control  the  maximum  temperature  T^,  This  clearly  will  be  reached  at 
y  =  0.  For  arbitrary  f(x),  the  only  general  way  to  find  the  x-coordinate  for 
the  maximum  temperature  is  to  calculate  and  plot 

/•a  oo 

Vo(x)  =  «(x> =  i  /  f<x)  +  2 Get)  cos  nr  tanh  rr 

J°  n=  1 

In  the  present  case  of  uniform  heat  generation,  it  is  apparent  that  Tm  will 
be  at  x  =  0,  under  the  center  of  the  emitter. 

If  we  define  <p  -  (T  -  Ts)/(q/k)  i  s  a  dimensionless  maximum  temper¬ 
ature,  then  the  desired  result  is  obtained  by  setting  x  =  0,  y  =  0  in  Eq,  (18) 
to  obtain 


I 


e.  .  nwx  . 

f(x)  cos - dx 

a 


(19) 


4(1) 


1 


(w/a)ir 


oo 

£4 

T  n 

n- 1 


sin  nir 


(r)  <“> 
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; 


Plots  of  calculations  based  on  this  result  are  shown  in  Figs.  3,  4,  and  5.  It 
should  be  noted  that  the  series  converges  slowly  for  small  w/a,  25  terms 
being  necessary  for  w/a  =  0.  1. 

If  the  emitter  occupies  the  whole  area  of  y  =  0,  which  of  course  is  not 
a  practical  device,  then  w/a  =  1  and  sin  mr(w/a)  =  0.  For  this  case, 
p  -  (l/2)(b/a),  which  is  the  solution  for  uniform  heat  flow  between  an  infinite 
parallel  source  and  sink.  We  call  this  p^(b/a).  For  b/a  >  1,  tanh  ir(b/a)  «s  1, 
and 


p  = 


(21) 


where 


1 

(w/a)ir 


sin  nir 


This  means  that,  for  a  large  device  thickness  b,  the  solution  is  the  uniform 
heat  flow  solution  plus  a  function  of  the  ratio  of  the  emitter  width  to  emitter 
spacing  w/a.  This  fact  is  illustrated  by  the  straight  line  portions  of  the 
curves  given  in  Fig.  3. 

A  parameter  of  interest  to  transistor  designers  in  the  past  has  been  the 
distance  at  which  a  neighboring  emitter  first  affects  the  thermal  behavior  of 
the  emitter  region  under  consideration.  If  we  consider  the  emitter  width  2w 
and  device  thickness  b  to  be  fixed,  then  the  variation  in  maximum  tempera¬ 
ture,  as  given  by  p,  is  a  function  of  emitter  spacing, as  shown  in  Fig.  5.  As 
a/b  decreases,  the  emitter  spacing  grows  smaller.  It  is  seen  that  the 
interference  between  emitters  has  effectively  vanished  for  a/b  >  2.  5.  This 
suggests  that  we  can  use  the  solution  given  by  Eq.  (18)  for  the  isolated 
emitter  solution  desired  in  Section  I,  if  we  only  set  a/b  >  2.  5.  The  value 
of  a/b  should  not  be  set  too  high  because,  as  mentioned  previously,  this 
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leads  to  slow  convergence  of  the  solution.  Consequently,  a  value  of  3  for 
a/b  was  chosen  as  a  compromise  to  give,  at  y  =  0  for  an  isolated  emitter, 


T  oo 

y=0  ~  a  =  1  +  3 

Z  .2(w/b)^n2 
n=  l 


sinimr(^) 


cos 


tanh 


nw 


(22) 


For  a  series  of  n  emitters  spaced  2a  apart,  we  must  add  the  solutions 
for  each,  substituting  x  ±  j2a  for  x  in  the  solution  for  each  particular 
emitter.  Here  j  is  the  integral  number  of  spacings  2a  by  which  each  partic¬ 
ular  emitter  is  to  the  left  or  right  (±)  of  the  emitter  under  consideration 
(that  is,  the  one  at  x  =  0).  If  we  perform  this  summation  and  solve  for  the 
temperature  at  the  center  x  =  0  of  the  emitter  under  consideration,  we  obtain 


*  = 


m 

T 


W-t 


7r2(w/b)  7 

1  n 


-j  cos  -j  nuj^j  tanh  (23) 


where 


j  =  0,  ±1,  •  •  • 

We  return  now  to  the  original  problem  of  an  infinite  number  of  emitters. 
In  designing  a  power  transistor  from  a  thermal  standpoint,  it  is  desired  to 
maximize  the  total  current  I  of  the  device  without  exceeding  a  maximum 
temperature  Tm  anywhere  within  the  device.  If  we  consider  a  device  of 
width  L,  then  the  total  number  of  emitters  present  is  L/2a.  Hence,  the 
current  in  any  single  emitter  is  i  =  I/(L/2a).  If  it  is  assumed  that  the  heat 
generated  beneath  any  emitter  is  q  =  iV,  then 

q  =  2IaV/L 
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q/k 


Fig.  3.  Maximum  Temperature  in  Power  Transistor  as  a  Function 
of  Ratio  of  Device  Thickness  to  Half  of  Emitter  Spacing; 
Heat  Assumed  to  be  Liberated  Uniformly  Under  Emitter 
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But,  from  the  definition  of  0, 


q  = 


k(T 


m 


T 

8 


therefore, 


I  = 


Lk(T  -  T 
m 

- ZaV0 


(24) 


The  maximum  total  device  current  results  when  l/a0  is  maximum. 

If,  for  example,  b  and  w/a  are  regarded  as  known,  then  this  will  be  at  the 
maximum  of  (b/a)/b0.  With  a  constant  \  defined  as 


Lk(T  -  T  ) 
m  s 

- 2F7 - 


then  Eq.  (24)  becomes 


I  -  v<b/a> 
0 


The  maximum  I  then  occurs  at 


81 

W5T 


90 

wiy 


=  0 


(25) 


or 


.  _  b/a  80  9  In  0 

1  =“T  WST  =  FTn  (E7a) 


(26) 
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In  other  words,  the  maximum  current  would  be  at  the  point  on  a  plot  of 
In  £  ve  In  (b/a)  where  the  slope  is  1.  As  might  be  expected,  no  such 
maximum  occurs  for  the  present  heat  distribution.  The  best  heat  dissipation 
will  be  obtained  when  w/a  =  1,  which  is  no  longer  a  practical  transistor. 
Having  a  minimum  possible  device  thickness  b  will  also  give  the  maximum 
heat  dissipation. 

The  foregoing  results  may  also  be  reached  in  a  more  intuitive  manner 
as  follows:  By  examination  of  Fig.  5,  it  is  seen  that  the  value  of 
(T  -  Ta)/(q/k)  increases  as  the  emitter  spacing  is  increased.  On  the  other 
hand,  more  emitters  can  be  placed  in  the  same  space  for  closer  emitter 
spacings.  Hence  the  maximum  heat  input  to  the  device  would  be  given  by  the 
maximum  of  q/a.  This  will  occur  at  the  minimum  of  $(a/b)  for  constant 
w/b.  If  this  operation  is  performed,  it  iB  again  found  that  the  optimum  occurs 
at  w  =  a,  for  which  the  emitters  are  all  touching. 

We  should  note  at  this  time  that,  in  the  preceding  analysis,  we  have 
assumed  Tfl  to  be  a  constant,  independent  of  the  total  device  heat  flux.  This 
is  not  true,  however,  because  the  device  itself  has  thermal  resistance  to 
the  flow  of  heat  to  the  surroundings,  and  T^  will  increase  with  the  total 
device  heat  flux  in  an  approximately  linear  fashion.  Nevertheless,  this  fact 
will  not  change  the  above  conclusions.  • 

For  non-uniform  heat  liberation  in  the  emitter,  however,  it  should  be 
possible  to  obtain  an  optimum  device  geometry.  In  any  case,  it  is  clear  that 
the  final  design  will  not  be  dictated  solely  by  thermal  considerations.  Elec¬ 
trical  performance  characteristics  and  available  fabrication  techniques  are 
prime  determining  factors  in  any  design.  The  effects  of  the  inevitable  design 
compromises  are  shown  as  the  squares  in  Figs.  3,  4,  and  5,  which  represent 
the  geometry  of  a  commercial  10-amp  power  transistor. 

C.  EMITTER  PARTLY  EFFECTIVE 

In  this  section,  we  consider  the  case  in  which  the  current  is  assumed 
to  flow  through  only  the  outside  section  wg  of  the  emitter,  but  uniformly  in 
that  portion.  This  condition  is  expressed  by 
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f) 


Fig.  5.  Maximum  Temperature  in  Power  Transistor  as  a  Function 
of  Ratio  of  Half  of  Emitter  Spacing  to  Device  Thickness; 
Heat  Assumed  to  be  Liberated  Uniformly  Under  Emitter 
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f(x)  =  for  w  -  <  x  <  w 


f(x)  =  0 


for 


0  <  x  <  w  -  w  and  w  <  x  <  a 
—  e  — 


Substitution  of  this  condition  into  Eq.  (2)  and  Eq.  (4)  yields 


Vs  ■  Tk-2  +  [ ""  ”(r)  •  “n  nTr(~x~^)] 

e  n=l 

(l)  8inh  nir(^T^)  8ech  n"(|) 


X  cos  nrrl  —  J  sinh  nit 


(27) 


We  will  leave  detailed  examination  of  this  result  to  the  reader. 

D.  OTHER  HEAT  DISTRIBUTIONS 

Solutions  for  other  current-heat  source  distributions  may  likewise  be 
easily  found  by  substitution  of  the  respective  functions  f(x)  into  Eq.  (2)  and 
Eq.  (4).  The  results  may  be  quite  complex,  but  the  principles  are  identical 
to  those  outlined  above  for  the  most  elementary  case. 

IV.  CONCLUSION 

The  internal  heat  transfer  problem  for  an  interdigitated  power- 
transistor  structure  has  been  solved  analytically.  The  relations  among 
current  distribution,  heat  generation,  and  temperature  distribution  have 
been  shown.  Application  of  the  results  to  the  most  elementary  case  illus¬ 
trates  usage  of  the  equations.  In  conjunction  with  considerations  of  other 
device  characteristics  and  current  fabrication  techniques,  these  heat  transfer 
results  should  permit  a  more  optimum  design  of  power  transistors. 
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NOMENCLATURE 


A 

a 

»n 

B 

b 

bn 

f(x) 

g(x) 

I 


j(x) 

Jo 

K 

k 

L 

m 

n 

q 


B 

T 


T 


m 


constant 

one-half  of  emitter  spacing,  cm 

nth  constant  in  solution,  given  in  Eq.  (4) 

constant 

thickness  of  collector  region,  cm 

Fourier  series  constant  defined  by  Eqs.  (6)  and  (8) 

1/k  of  heat  liberation  rate  at  any  point  x  under  emitter,  °C/cm 
temperature  (°C)  at  any  point  x  on  y  -  0 

total  current  of  all  emitters  per  unit  length  of  single  emitter, 
amp/cm 

total  current  for  one  emitter  per  unit  length  of  emitter, 
amp/cm 

2 

current  density  at  any  point  on  emitter,  amp/cm 

constant  in  Eq.  (14),  amp/cm^ 

constant  in  Eq.  (11) 

thermal  conductivity,  watt/cm  °C 

length  of  device,  cm 

total  number  of  emitters  (Eq.  23) 

index,  0,  1 ,  2,  ... 

heat  liberated  at  each  emitter  per  unit  length  of  emitter, 
watt/  cm 

integral  constant,  0,  1,  2,  ... 

temperature  at  (x,  y),  °C 

maximum  temperature  in  transistor,  °C 
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NOMENCLATURE  (Continued) 


T 

V 


a 


w 

w 


e 


x 

x 


e 


y 

y 


* 


w/a 

b/a 

sinh  x 


cosh  x 
sech  x 


temperature  of  heat  sink  at  collector,  °C 

effective  potential  drop  for  current  flow  in  volts,  so  that 
q  =  iV 

one  half  of  emitter  width,  cm 

one  half  of  effective  emitter  width,  from  which  current  flows, 
cm 

distance  from  center  plane  of  emitter,  cm  (see  Fig.  2) 
distance  from  outside  edge  of  emitter,  cm  (w  -  x) 
distance  from  top  surface,  cm  (see  Fig.  2) 


Lk(T  -  T  ) 
'  m  s' 

- 2F7 - 


T 

8 


q/k 


ratio  of  device  width  to  emitter  spacing 


oo 


(w/a) 


a)TT  n 

'  n=l 


ratio  of  emitter  width  to  emitter  spacing 
twice  ratio  of  device  width  to  emitter  spacing 


x  -x 

e  -  e 


x  -x 
e  +  e 


1  /cosh  x 
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NOMENCLATURE  (Continued) 


cosech  x 
tanh  x 


coth  x 


1 / sinh  x 

sinh  x/cosh  x 

cosh  x/ sinh  x  =  1/tanhx 
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